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Abstract. A new set of nonlinear coupled equations is derived in the context of small 
amplitude limit of the general wave equations in a fluid type warm electrons/cold ions 
plasma irradiated by a continuous laser beam. This limit is proved to be integrable 
by means of the spectral transform theory with singular dispersion relation. An exact 
asymptotic solution is obtained. This model accounts for a nonlinear mode coupling of 
the electrostatic waves with the ion sound wave, and is shown to be unstable and does 
not propagate any stable small amplitude solution. This instability is understood as a 
continuous secular transfer of energy from the electrostatic wave to the ion sound wave 
through the ponderomotive force. The exact mechanism of this transfer is exposed. 
The dynamics of this energy transfer results in a singular asymptotic behavior of the 
ion sound wave which explains the low penetration of the incident laser beam. 


PACS numbers: 02.30.1k , 52.35.Fp , 52.38.-r 

Keywords : inverse scattering transform, singular dispersion relation, nonlinear couples 
waves, laser-plasma interaction 


Submitted to: Nonlinearity 


New example of integrable nonlinear coupled equations... 2 

1. Introduction 

Recent progress in laser generation has aroused a great interest in the problem 
of laser-plasma interaction, which is the source of various nonlinear effects such as self- 
focusing, self-modulation, hlamentation and several parametric instabilities [U El [3l HI |5] . 
These phenomena have applications in many advanced physical events such as inertial 
conhnement fusion, higher order harmonic generation. X-rays sources and particle ac¬ 
celeration pmiHiE]. 

This work regards the general theory of nonlinear coupled waves in plasmas, in the 
context of a two-components non magnetized fluid-type plasma. In a plasma produced 
and irradiated by laser beam, the electrostatic wave (EW) induces scattered waves by 
means of different types of stimulated emission mm- We are interested here in the 
process of stimulated Brillouin (back) scattering (SBS) for longitudinal EW. The result¬ 
ing incident and reflected high frequency EW can drive some low frequency waves of 
acoustic type, such as the Ion Sound Wave (ISW), by means of the ponderomotive effect 
on the local charge density. In the present paper, it is shown here that this nonlinear 
effect is at the origin of plasma instability. 

It is well known that the model of nonlinear propagation of ISW, when the pon¬ 
deromotive effects are neglected (no mode-coupling), is the famous Kortoveg de Vries 
(KdV) equation [TUI ITT] . This equation has the advantage of being integrable, which 
allows to describe the solution on the basis of its ’’nonlinear Fourier spectrum” (see f.i. 
[T2] and m) . When the ISW is decomposed into his harmonic modes, the KdV equa¬ 
tion maps into the nonlinear Schrodinger equation (NLS), an other integrable model 
[iniiiiiiniiia]. Both models propagate soliton solutions which provide the general fluid 
equations with stable small amplitude localized coherent structures. But these models 
ignore the possibility of energy transfer from the long wavelength domain to the dissi¬ 
pative short wavelength range. 

This energy transfer is considered by including in the momentum transfer equation 
for the electrons, the poderomotive force (low frequency effect of the high frequency elec¬ 
trostatic held) which results on the nonlinear mode coupling. In this approach Karpman 
m derived an integrable model, which however does not account for dispersion or non¬ 
linear dynamical effects (ISW nonlinearity). Still this model has proved to be very useful 
and has been slightly expanded in [15] and [16]. All these models and their stability 
properties are discussed in mi- 

Later, M. Casanova and ah ng, have shown that the wave dynamics in plasma 
is well described by a model where the ISW is nonlinear and is coupled to the elec¬ 
trostatic held through the poderomotive term. In particular, the ISW nonlinearity has 
been shown to play a crucial role in the Brillouin rehexivity. 
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In this paper, starting from the most general one-dimensional fluid equation for 
a plasma with a Maxwellian velocity distribution [10], a model for the nonlinear 
propagation of the ISW coupled to the EW (the pump wave and the wave reflected 
through SBS) is established. This model is obtained as an exact asymptotic limit in 
the regime of ’’warm electrons-cold ions” and small amplitude of the pump wave. The 
technique used is the well known multisclae expansion naEi], and the model exposed 
in this paper, appears then as universal in the sens given in 121 .This model is a new 
nonlinear evolution equation with unexpected properties. The integrability permits to 
obtain general theorems about the nature and the behavior of the solution. In particular, 
it is shown that the EW transfers its energy to the ISW through a highly nonlinear and 
unstable process : any inhnitely small initial disturbance of the ISW asymptotically 
evolves towards a singular (soliton-like) solution. 


2. Derivation 

Let consider a linearly polarized electrostatic held vector E = zE{z,t) with a 
frequency close to plasma frequency ojq = dvre^no/me with some broadening represented 
by the following wave paquet 

/ + 00 

dujE{u}, z,T)e~^‘^^ . ( 2 . 1 ) 

■CXD 

Note that mo and —e are respectively the mass and the charge of electrons and no is 
the ambient ion density. E {u,z,t) is the slowly varying envelop of E{z,f) where r is 
the slow scaled time dehned below [2T^ 


The low frequency effect on the electrons due to the high frequency held E, results 
in the ponderomotive force obtained by considering the real part of [Sz{dE/dz)] around 
the average position < > where 5.^ is the complex valued displacement and solution 

of 


me{Sz)tt = eE . 


( 2 . 2 ) 


The ponderomotive force reads then 


'‘+00 


fp = 


2mp 


dz 


du u 


-2 


E{u,z,t) 


(2.3) 


Using the ratio of electron and ion masses ei = {me/miY and the electron Debeye 
wavelength A|, = KsTe/ATmoe"^ (with the electron temperature and Kb the Boltzman 
constant), a dimensionless set of variables can be set as follows : 


z' = \bZ , T = eiOJot. 


(2.4) 


in these variables, we have the usual system of huid-type equations [TO] with the only 
additional term fp in the momentum transfer equation for the electrons |22j . 
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Moreover, the Maxwell equation for the EW is obtained on the basis of the 
dispersion relation [10] 

u^ = ujI + (2.5) 

where Vre = ^d^o, is the thermal electron velocity, k is the wave number and ujp stands 
for the plasma frequency 

+ q,) . (2.6) 

Qf. is the fractional change in the electron density of average value uq, namely 

Ue = no [1 + ge(^,f)] • (2.7) 

Now, redehning the helds as 


& = E 




mi 


1/2 


(2.8) 


where (j) is the electrostatic potential and Vi the ion velocity, the momentum transfer 
equation for electrons reads [iniES] 

2 

= 0 , 


d(j)' 1 dqe d 

dz' (1 + ge) dz' dz' 

and for the ions 


p+oo 

/ du 

E\u},z',t) 

J —oo 



dv', . dv'y d(j)' 

—- + V —- = - — . 

dt' * dz' dz' 

The above system has to be completed with the continuity equation 

^ + ^[(i+9<)<'a = o, 

and the Poisson equation 
d^f)' 


dz'^ 


= , 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


where g* in 12.111 and l2.12l is the fractional change in the ion density n* = ne(l + g*). 

Let dehne a slowly varying set of variables in the comoving frame of the ISW at 
the sound speed Cs = {KBTe/miY^'^ = eiUoXn ; 

f = ei{z'-t') , r=(ei)^'^^f. (2.13) 

Now, by expanding the electrostatic held 0, the density variations ge and g* and the 
scaled ion velocity n', in powers of e as 

t>' = + ePf> + 0(ej) (2.14) 

and by dehning S as follows 

& = {eif^ [S + 0{e,)] , (2.15) 

the set of equations 12.9112.12] after integration of 12.91 gives at the hrst order 


(2.16) 
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and at the second order, by eliminating and , one obtains the following 

evolntion for r (dehned in l2.16p : 


2 r^ + 2rr^ + 


'• + 00 


-d, 


\£\ 


(2.17) 


On the other hand, the Maxwell eqnation from 12.51 reads 


dt'^ 



E{z,t) ^ -Wg [1 + q^a.i)] E{z,t) . 


(2.18) 


Using the change of variable 12.41 the scalings 12.131 12.151 and the expansion 12.141 the 
eqnation 12.181 gives at hrst order 


~ q)^ — 0 5 


(2.19) 


where the parameter u is dehned throngh the relative difference between the sqnared 
freqnency u of the applied laser beam and the sqnared plasma freqnency Uq : 

= 3 . (2.20) 

^0 

The system 12.17112.19] representing the nonlinear propagation of the EW conpled to the 
ISW does not acconnt for SBS and has been shown to be not integrable [ 22 ]. To take 
SBS into acconnt, one has to allow the electric held envelop E to contain simnltaneonsly 
at hrst order, a left-going wave exp [i{uJiT + ki)] and a right-going wave exp [i{u} 2 T + ^ 2 )] 
with small slowly varying amplitndes 6201 (.^,r) and 6202 (.^,r) respectively, 62 being a 
small parameter snbmitted to the constraint 12.291 The wave nnmbers modnli are closed 
to a same average valne ke, namely : 


ki = ke + 0 ( 62 ) , k 2 = - ke + 0 ( 62 ) 


( 2 . 21 ) 


which means 


E{f,t) = 62 




i(a;ir -|- kiz) 


-|- 02 e 


i(a;2r k2z) 


( 2 . 22 ) 


and the freqnencies are related throngh (resonant scattering conditions) 


uji — LjJ2 + ojg , ki — k2 kg (2.23) 

where Ug and kg stand for the ISW parameters. In the frame {f, r), the dispersion rela¬ 
tion reads Ug = ^{kg)^ . 


In the laboratory frame {X,T) with the scaled variables [19] dehned as follows 


X = e 2 


^ T 


T = 3 ( 62 )^ kg T , 


(2.24) 


we consider r{^, r) (dehned in I2.16P to be a snperposition of its harmonic modes 
expanded in powers of 62 as the formal series 

N 

r = {€ 2 ? qo + Yl ^" 2 )" qn + ^20 + c.C. (2.25) 

n=l 




























New example of integrable nonlinear coupled equations.. 


with 


M 


= > g.,(X,T)( 62 r + 0(6"+^) 


j 


n = 0,N . 


(2.26) 


For N = 2 and M ^ 3, by means of the standard mnltiscale techniqnes |19l 1201 122] 
i.e. by inserting 12.221 and 12.251 into 12.171 and I2.19[ then examining the relevant orders 
in the expansions in powers of e 2 and averaging over the fast oscillation exp(ia;ir) 
(which means : look at the coefficients of exp(ia;sr) only, known as the ’’rotationg wave 
approximation”), one obtains for 

9=X«o(V,T) , (2.27) 


3iks 

the following system of conpled eqnations 


2/cp 


+ 2^xx - (1 y dC,aia2 


(2.28) 


0-1,X N iC®i = QO 2 
02,x ~ K 02 — QOi 

the left going wave of amplitnde ai(X, T) resnlts from the incident (pnmp) laser beam 
with wave nnmber ( and initial given amplitnde, say 7l((C). The right going wave of 
amplitnde a 2 (X, T) is the wave reflected by SBS. 


Note that the two mnltiscale analysis with the small parameters ei and 62 applied 
in this section have to be consistent, which is gnaranteed by the ordering 


ei < 62 < 62 -C 1 . 

To the system [2.281 we adjoint the following initial/bonndary valnes 
(q(X,0)€L\ 


lim ai{(,X,T) = A{C)e-‘<^ 

X^+OO 

lim a 2 {C, X,T) = 0 

X —¥ — 00 


(2.29) 


(2.30) 


Thns, the system of conpled eqnations 12.281 completed by the initial/bondary conditions 
12.301 is the central point of discnssion in the following sections. 


3. Method of resolution 


The system [2. 281 with the initial/bonndary conditions 12.301 belongs to a general class of 
integrable nonlinear evolntion [23] with a singnlar dispersion law |23l |2Tl |28l |29]. This 
system is a special case of the following system for the 2 x 2 matrices fi{k,x,t) and 
Q{x,t) : 


Qt - ^osQxx + = i [0-3, (/iCTs/i ^)] 

/Ta; + ik [0-3,/x] = Qp- 


(3.1) 
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where 


and 


Q = 


0 q 
q 0 


= 


/^ll P'12 
p21 P22 


0-3 = 


1 0 
0 -1 


(/iCTa/i ) ^ ^ JJ dX A dX g{X,t)p{X,x,t) as p ^{X,x,t) , 

for any given distribution g{k,f) in the complex plane. We have also 
dX A dX = —2i dXR dXj for A = + iAj 


(3.2) 

(3.3) 

(3.4) 


The above system has been shown to be integrable [211 [28l [29]. It is related to the 
general theory of the spectral transform [El Ell in the context of the singular dispersion 
relation [ 22 ] • A complete proof and details are given in the Appendix. 


One can summarize the method of resolution by the following set of linear steps : 

Step 1 : Given q{x,t)) on L^(M), solve the following Volterra equations for the 
vector /i+((C, X, 0) and /i“((C,X, 0) with components pi and p 2 


< 


/•+°o 



n+(X> 

pX 

= 1 - 

/ dX'qpf 

p 7 

= 1 - 

/ dX'qpf 

A 


Jx 



Jx 

\ 


fX ’ ^ 



r+00 

A 

= 1 - 

/ dX'qvf 

^2 

= 1 - 

/ dX'quf 

< 


J —00 

K 


Jx 


with u = /iexp(2i^X) . 

The vector pN (respectively p~) is homomorphic in A((C) > 0 (respetively in 
A(C) < 0) and the discontinuity p~^ — p~, as C crosses the real axis leads to the following 
Riemann-Hilbert problem for C ^ ® 

P~^{C + iO) - p~{C - iO) = R{C + iO) e^“’^ai/i“(C - iO) . (3.6) 

The function R{( + i0) is called the nonlinear Fourier Transform of q{X, 0) and ai is the 
usual Pauli matrix. R is computed out of the Volterra equation 13.51 for given q{X, 0) by 

/ +00 

dX'g(X',0)/i+(C,X',0)e-2*^^'. (3.7) 

■00 

This hrst step solves the direct problem, namely the determination of R{C) from 
q{X,0). 

Step 2 : When q{X,T) obeys the nonlinear evolution equation of 12.281 the time 
dependance of R is given by 

R(C,T) = R(C)e-^^ (3.8) 

with 


(3.9) 
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where the slashed integral denotes the Cauchy principal value. It worth remarking that 
6 has a positive imaginary part 6i = (7r/12)a;o l^(C)l^) which implies thes exponential 
growth of R{( + iO,T) in time. We shall see that this fact is the very reason of the 
instability for the system 12.281 

Step 3 : Solve the Riemann-Hilbert problem 13.61 in which R{(, T) is given in terms 
of i?(C,0) by means of 13.71 The solution is obtained from the following Cauchy-Green 
integral equation (omitting the (X, T)-dependence) 

where C represents any contour in the upper half plane (7y((C) > 0) equivalent to the 
real axis. Equation 13.101 gives the solution while the solution fi~ X,T) is 

obtained out of EH 

Step 4 : Finally the solution {q, 01 , 02 ) of the system [2.281 is obtained from the 
function /i through 

f«(Vr) = 2ifi«(V.T) 

\ a,(X, T) = .4^+(C, X, , i = 1, 2 

where fb^^'>{X,T) denotes the coefficient of in the Laurent expansion of /22(C) 

This last step achieves the resolution of the initial/boundary value problem for the 
system 12.28112.301 

Note that the equations 13.51 and 13.61 imply 

lim = (3.12) 

X —^-j-oo 

Hence, R(C + iO) is in fact the recflection coefficient of the potential q at any hxed time 
t. It is easy to verify that the transmission coefficient T, dehned by 

lim /i+ = r, (3.13) 

x—^—00 

and the reflection coefficient R obey the unitarity relation 

R(C + iO) R(C - iO) + r(C + iO) r(C - iO) = 1 . (3.14) 

4. Conservation laws 

The hrst usual consequence of the integrability of a system is the existence of an 
inhnite sequence of local conservation laws. 
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The explicite calculation of conservation laws requires the Lax pair Q{k) and W{k) 
Here, is the (singular) dispersion relation (this choice is justihed in the Appendix. 
See 16.351 and 16.361) 


kl = ifc^cTs H—^ 


dX A d\ 


2i7r JJq X — k 


g{X)a3 


and W can be computed with the methode exposed in [22] 


2 I ; I ^ ^ 


W{k) = -iask^ + Q (Qx-Q )-7^ 


dX A dX 


2i7r .//j- X — k 


(4.1) 


g{X)'ipa3i) (4.2) 


0 q 


where Q = ( I ^1 and tb is dehned in 16.21 and 16.31 
\q OJ 

Using the compatibility condition 
Ut-Wx + [U ,W] = 0, 


(4.3) 


the inhnite sequence of conservation laws can be obtained by equating the powers of 
{2ik)~^. The hrst law of this sequence reads 


drlql =d_ 


'X 


^ {qxq - qxq) + ^ ^ 


with 


'‘ + CXD 


J = 


dC [®i(C T iO)®i(C ~ iO) T ®2(C T iO)®2(C ~ iO)] • 


(4.4) 


(4.5) 


Now by integrating the equation 14.41 with respect to X using I2.3UI 13.61 and 13.141 one 
obtains 


A 

&T 


'•+00 


'‘ + CXD 


dX' \qf = dCl^lTR(C + iO)i?(C-iO) 


(4.6) 


where |A'| = (a;o/12) |A|. 


Due to 13.81 the product R{( + iO)R{C — iO) is time-independent and therefore, the 

/ -l-oo 

dX' |g|^ transferred from the EW to the ISW, grows linearly in time. This 

-OO 

causes the instability of the ISW which is going to be described in the next section by 
evaluating the time-asymptotic solution of q{X,T). 


5. The time asymptotic solution 

As show in equation 13.81 the reflection coefficient R{( + iO,T) has an exponential 
growth which is characteristic of integral evolutions having a singular dispersion law 
[23l 1211125] . Following [25], the asymptotic behavior as T ^ -|-cxd, of the solution of l3.10l 
is evaluated by the saddle point method. 
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Let consider an arbitrary initial local distribntion q{X, 0) on the charge density by 
considering some generic reflection coefficient 

ro 


R{C + iO , 0) — 


(5.1) 


(C + iO)" ’ 

where tq is a real constant and n > 1. Let consider also a very sharp (i.e. 
monochromatic) inpnt laser beam, so that the evolntion 13.81 is rednced to 

(5.2) 

'* + 00 


i?(A) ~ tqA ”exp iT(Y — 2A^) 
L A 


for Q=(A) > 0 and 7 = / dC,' |A'((C')|^. From 15.21 the saddle point as T ^ +cx), is 


_ 1 + i\/3 7 (1/3) 

+ -^—( 4 ) 


(5.3) 


One may choose for the contonr C in 13.101 any contonr eqnivalent to the straight line 
parallel to the real A-axis passing thonght this saddle point. This contonr, eqnivalent to 
the real axis, is precisely the critical path. Then the asymptotic behavior of the solntion 
of 13.101 satisfies the following algebric system (9'(C) > 0) 

1/12(0 y' ■ 

where the constant B is 


1\ 2^-(i/2) 


2i\,x + ir (7 - 2Xi) 

A.o 




B = 


71 ro 


6 (A,)-" 2i7r • 

The solntion of the above algebraic system gives the behavior of and 

fi'^{C,X,T) throngh [3.61 and hence, Oi and 02 throngh 13.111 Namely 




A.' i-|T(A:,r)|‘ 




(5.6) 


A.' i-|T(A:,r)|- 


with 


where. 


T(X,T) = 


A. 


As — As 2i7r [Ta;"(As)_ 


-27r 


a(As) e 


i<F 


$ = 2iAsX + 2i(A2 - ^)T 

and u” stands for the second derivative of a;(A) = 2i 

A. 


AA' 7 


(5.7) 

(5.8) 

with respect to 


Finally the time asymptotic behavior of q is obtained throngh 13.Ill 


lim 

T—>.+00 


q{X,T) 


-2^A(As) 


ei0 

sinh p ’ 


(5.9) 
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0 = 23fi(A,) [X + 63fJ(A,)T] + 0o 


(5.10) 


and 



where the constants a and 0o are related via 



(5.12) 


The consistency of the solntion 15.61 and 15.91 can easily be verified by replacing these 


solntions into the initial system 12.281 satisfying the initial/bonndary condition 12.301 


6. Conclusions and discussion 

In conclusion, we make the following comments : 

(i) Although the asymptotic behavior of q{X,T) 15.91 is a singular expression, q{X,T) 
itself, is not singular. Thus the direct consequence of 15.91 is the growing amplitude 
of q{X,T) which is the signature of the energy transfer. 

(ii) The singular point in 15. 9 1 travels with the asymptotic velocity —6 3fJ(As). Hence, the 
solution q{X,T) accumulates energy in the region where the laser beam is applied, 
which accounts for the low penetration of the laser in the plasma. 

(hi) The basic system 12.281 is obtained as the small amplitude limit of the fluid type 
behavior of plasma and the Maxwell equations. Therefore, it is proved here, that 
this system does not propagate stable small solutions when there is a balance 
between SBS emission and the ISW nonlinearity. It worth remarking that, by 
neglecting the ponderomotive effect (which means no right hand side in the 
evolution equation in I2.28p . one would obtain an equation for the ISW where any 
initial disturbance would eventually disperses away as T —)• cx). Therefore, while the 
nonlinearity in the left hand side of the evolution equation in l2.28l is a mechanism for 
the saturation of the (linear) instability of the ISW, the nonlinear ponderomotive 
effect of the electric held is the mechanism responsible for the plasma instability. 

(iv) It has been admitted that q{X,T) G L^(M) which implies in particular that 
\q{X, T)| —>■ 0 as |X| —)■ oo. It is interesting to wonder whether different asymptotic 
behaviors might change the stability properties of the system. For instance, the 
large thermal conductivity of the electrons might act as a saturation mechanism. 

(v) In many physical situations the length of the plasma is hnite. Then, arises the 
question of the properties of the system 12.281 on a hnite X-interval. It is possible, 
as in the study [18], to observe in this case a chaotic behavior. 
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(vi) It is well known since the pioneering work of Zakharov [2S], that the ronte 
to Langmuir turbulence passes through the ’’collapsing of the EW”. i.e. the 
formation of local singularities of the wave amplitude (and hence also of the ISW 
amplitude named cavitons). This result has been originally obtained in 3 dimension 
with spherical symmetry and later, in one dimension 123 but with higher order 
nonlinearity having no physical origin. Here we have a model equation, derived 
from the basic hydrodynamic-Poisson-Maxwell system equations, whose solution 
generically evolves through a singular solution. 

(vii) Note hnally that the monochromatic approximation led to the evolution 15.21 and 
subsequently to only one saddle point and one singularity in the time-asymptotic 
behavior of the ISW. Releasing the monochromatic approximation, will introduce 
more and more saddle points and hence, more and more singular points or 
asymptotically collapsing waves. 


( 6 . 1 ) 




U — —k Q 


Q = 


( 6 . 2 ) 


Appendix : Proof of the method of resolution 

The differential equations 

ai,x + iCoi = qa2 
0‘2,X ~ iC®2 ~ O'®! 

extracted from the system 12.281 is the vectoriel form of the Zakharf-Shabat spectral 
problem [301 El] that can be written in the general matrix form [23 

0 q 
r 0 

The case 16.11 is recovered by the reduction r = q, which is going to be adopted from 
now on. Here -0 is a 2 x 2 matrix, built with two independent column vector solutions 
("01 ,02) and is completely determined by giving its asymptotic behavior. Then the set 
of the differential equations 0a; = f/ 0 can be equivalently written as a set of Volterra 
integral equations. For convenience, let write these equations for the matrix qi{k,x), 
dehned by 

fi = 0exp(i/c(T3x) ^ /j,x = ik[n , as] -f Qp.. 

Two dependant solutions and p~ can be dehned through 

/ill ~ 1 


(6.3) 


/ill ~ 1 

At = -j 

1^2 ~ 1 ~ 


—oo 

P+CXD 


' dx' q P 21 

X 

d0g/i+ 

dx' q /i+ 


'‘ + 00 


'» + CxD 


'‘+00 


dx' q /i ^2 


1^21 — / 

J X 

1^2 = -/ 

1^22 ~ 1 ~ 


dx' q P 21 

dx' q /i 02 
dx' q /i 0 


(6.4) 


— OO 

pH-00 
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One can easily verify, by using the Leibnitz formulas, that /x satishes [6.3l In the reduction 
r = q, one has 

Q = aiQai , 
and it is easy to prove that 

p.'^{k,x) = (Ti p.~{k,x) (Ti . (6.6) 

The set of Volterra equations 16.41 can be mapped into a Riemann-Hilbert problem as 
follows : 


/i’*' (respectively /i ) is holomorphic in the complex upper half plane ^{k) > 0 
(respectively in the complex lower half plane Q'(fc) < 0). Let write 

f u'*' in 'N(k) > 0 

ia = < _ . (6.7) 

[ /X in ‘N{k) < 0 

The matrix /x is holomorphic all over the complex plane except on the real axis where 
it stands a discontinuity. Writing the Riemann-Hilbert problem consists on expressing 
this discontinuity in terms of /x^ (or /x“). 


To do so, let Compute 

D{k, x) = [frf{k + iO , x) — frf{k — iO , x)\ , xgM (6.8) 

out of the integral equations 16.41 and obtain the following Volterra equations for Di and 
D 2 , components of D. 



where 



dx' q{x') fiii{k,x') 


%k) = 0 + . 


( 6 . 10 ) 


An integral equation having the same Green function as 16.91 can be obtained readily out 
of I6.4[ Indeed, the vector 

L = frif {k + iO , x) a'^(k) , A; G M (6-11) 


is also a solution of 16.91 General theorems about integral equations [32] allow to prove 
that 16.91 has only the trivial solution a = 0. Then we have D = L, which gives the 
following Riemann-Hilbert problem 

iaf{k -f iO , x) — iaf{k — iO , x) = a^{k) + iO , x) . (6-12) 

A similar calculation with 

x) = [/X 2 (fc-|-iO , x) —/x^(fc — iO , x)] , xgM (6.13) 
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leads to the second Riemann-Hilbert problem 

/iJ (A; + iO , x) — iaf{k — iO , x) = a~{k) nf{k — iO , x) . (6.14) 

Using 16.41 we have 

/ +CXD 

dx' q{x') fif 2 {k,x') = a~^{k) , = 0“ (6.15) 

— OO 

Then 16.121 and 16.141 can be written in the following matrix form 

- fi-= {fif , fif) S . (6.16) 

with 

, keR. (6.17) 

Note that , as /i in 16.61 the matrix S satishes the reduction relation 

S{k,x) = —(Ti S{k,x) (Ti . (6.18) 

The function a^{k) is called the reflection coefficient. a~^{k) and q{x) are equivalent 
in the sens that q{x) being given, one solves the Volterra equations 16.41 and computes 
a~^{k) through 

a~^{k) = lim fi 2 i{k,x)e~‘^^’^^ , (6.19) 

a:—>-+oo 

which solves the direct problem. 

The inverse problem consists on constructing Q{x) from a given S{k,x) and solve 
the Riemann-Hilbert problem l6.16l to get p.{k,x). Once /x is obtained, then q{x) is com¬ 
puted in the following way : 

Write the Laurent expansion 

{ - 1 - (1) 1 - (2) 

/Xii = l + ^/iii + +•••• 

- _ 1 - (1) , 1 - (2) , 

7^21 - j^k‘21 + ^ 2/^21 + •••• 

obtained from the Volterra equation 16.41 by integrating by parts 
16.31 and use the Liouville theorem to obtain 

q = 2i /i 2 i 

which gives q{x) out of fi{k,x). 

To solve the Riemann-Hilbert problem, it is convenient to use the “DEAR problem” 
formulation : 

^ ^ i [m V(x,) - r(X,)] (6.22) 

where 

d V ^ ^ i 

dk 2 dkji dkj ’ 


( 6 . 20 ) 

. Then, insert 16.201 in 

( 6 . 21 ) 


S{k, x) = e 


_ —ifcc73S 


0 

a~^{k) 


—a+{k) 

0 


ikcr^x 


(6.23) 
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k = kji + ikj , h"*" = — iO) , S =5(/c7- + i0), 

6 being the usual Dirac distribution. 


(6.24) 


Thus the equation 16.161 becomes 
dp, 


^ = uK , U = - \ 

dk ^ 2 V 0 

And the generarized Cauchy forrmula reads 

1 




V 


dX A d\ 
X-k 


(6.25) 


/i(A, k)R{X, k) (6.26) 


with T) the complex plane. 


The Laurent expansiou l6.20l allows the computation of the hrst integral iu l6.26l and. 
taking into account the analytic properties of p, the second integral of 16.261 is reduced 
to an integration over the real line. Consequently, for the hrst column vector pi with 
the reduction of 16.61 we have 

( 1 \ 1 r+oo i\ _ 

oj+—y -^a,pfCKx)a+{X)e^^^N (6.27) 

A similar equation holds for the second column vector. 


Hence the inverse problem (construction of q{x) from a~^{k)) is resolved by solving 
the Cauchy-Green integral equation and by calculating q{x) out of 16.211 

Notice that the preceding formalism remains valid if q{x) is assumed to depend also 
on a real external parameter t (time). Then, the eigenfunction p{k,x) and the spectral 
transform R{k, x) will also depend on t. 


On the next step, we must construct integral evolutions out of a simple choice of a 
given time dependence of the spectral transform R{k,x,f). 


We start with the following generic DEAR problem 
^ /i(fc) = p{k) R{k) 

p{k) = l + 0{^) , |/c| CO 
and ask for R an {x, f)-dependance through 

Rt = [R, 

Rx = [R, A] 

where A and D are given distributions of A; G C, functions of x and t 


(6.28) 


(6.29) 
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It is a simple task to check the following relations 
d dA 

^ ( -1 o -1^ -1 

By integrating the above equations, one obtains 

d\ A dX 


ia^{k,x,t) = 


p,t{k,x,t) = 


1 ff dX A dX ,dA(X,x,t) 1 ,^ , 

U- — II — - —ia{X,xN )— iX,x,t) 


2711JJ c X-k 


dX 


X 


/i(fc, X, t) + p,{k, x, t)A{k, x, t) 


1 ff dX A dX dXl(X,x,t) 1 ,, , 

V- — II - - -fi{x,x,t )— ’ V 


27ri 


c 


X-k 


dX 


X 


(6.30) 


(6.31) 


p.(k, x, t) + fi{k, X, t)fl{k, X, t) 

in which the matrices U and V (the ’’constants” of integration of the d operator) are 
given by 

U{k,x,t) = —Polk [ia{k, x,t)A{k, x,t)ja~^{k, x,t)~\ 

V{k,x,t) = —Polk [fr{k, x,t)il{k, x,t)fi~^{k, x,t)j 
where Polk [....] means ’’the polynomial part in k of 


(6.32) 


To obtain the equation 16.31 we choose and adopt from now on 
A = ifc (Ts . 

The compatibility condition = k-tx when dA/dk = 0, leads to 
Ut{k)-V,{k) + [U{k),V{k)] = 

dX A dX TTfw —i/\n 

^ ^ U{X)-U{k) , fi{X)—^fi (A) 


27ri 


Now, with the choice 


f2 = ik^as H—^ 


dX A dX 


2i7r JJq X-k 
the time evolution 16.341 is transformed into 




Qt “1“ ^ 

which is the evolution equation in 16.311 


(6.33) 


0-3 , ^ // dX AdXg{X)ia{X)a3ia ^(A) 


(6.34) 

(6.35) 
, (6.36) 


Therefore, it has been proved that when Q{x,t) evolves according to a nonlinear 
evolution, its spectral transform R{k,x,t) related to a^{k,t) by 16.171 evolves in time 
according to 16.301 with Q{k) given in 16.351 


In the reduced case r = q we have the relation 16.181 therefore the explicite time 
evolution of the reflected coefficient a = a~^{k) is given by 

dXjidXj 


(^t{kR,t) = 2i 


K-- 


TT i/c A - {kR - iO) 


TTPdiX) 


a{kR,t) . 


(6.37) 
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The example of the physical moded studied in the section 2 justihes the most intresting 
following subcase : 


ITT 


^?(A) = -h(A,)h(AH-MI^(MI 

0 

Indeed, in the case 16.371 gives after integration 


a{k, t) = a(fc, 0) exp 


2i ( k 


'• + CXD 


+ 6 


dk 




(6.38) 


(6.39) 


, ko-k J 

Clearly the function a{k) has an essential singularity at k = k^. Hence a{k) is not 
dehned for k = k^ unless 0 (^ 0 , 0) = 0, which actually required by the system 12.281 
Indeed, for g(x, 0) vanishing at both ends of the x-axis, the consistency requires that 
01(22 vanishes too. Therefore, as a consequence of equations 16.141 and 13.111 a{k) must 
vanish for k = k^. This condition actually determines the parameter ko from the initial 
datum g(a;, 0), which is the solution of the equation 

r+oo _ 

Q((fco,0)= / dx'q{x') pL 72 {kQ^x\t)) =0. (6.40) 


Hence, it is a property of the system 12.281 that the initial value q{x, 0) of the ISW 
determines the small correction of the wave number of the scattered ESW. 
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